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Abstract 

The spectral properties of the 1-D Hubbard model are obtained from quan- 
tum Monte Carlo simulations using the maximum entropy method. The one- 
particle excitations are characterized by dispersive cosine-like bands. Veloc- 
ities for spin- and charge excitations are obtained that lead to a conformal 
charge c = 0.98 ± 0.05 for the largest system simulated (N = 84). An exact 
sum-rule for the spin-excitations is fulfilled accurately with deviations of at 

most 10% only around 2kp . 
PACS numbers: 75.10.Jm 



1 



Typeset using REVTgX 



Since Anderson [1] proposed that the high-T c superconductors (HTS) should be con- 
sidered as Luttinger liquids [2], a great deal of interest was focused on such systems. The 
best known Luttinger liquids are the Luttinger and the Hubbard model (HM), both in one 
dimension (1-D). In spite of the fact that both models can be solved exactly [3,4], very little 
was known about their dynamical properties until recently. 

In the case of the Luttinger model, the spectra for collective excitations can be easily ob- 
tained by bosonization [3], whereas the one-particle properties were calculated only recently 
[5,6]. However, since lattice effects are neglected in this model, the results obtained give 
only the asymptotic behavior for vanishing excitation energy. Being this region the most 
difficult to be accessed by experiments, further progress is necessary in order to clarify the 
situation for the HTS. 

Lattice effects are contained in the I-D HM, that has an exact solution by Bethe-Ansatz 
(BA) [4]. A number of authors succeded in extracting from BA information about spectral 
properties [7-11], however, many of these results are limited to special situations (e.g. half- 
filling, U — > oo and/or one-hole doping, or u = 0). Only recent progress achieved in the 
frame of conformal field-theory [12], led to the asymptotic properties of correlation functions 
irrespective of the coupling constant and doping. Complementary to these achievements, a 
Landau-Luttinger liquid theory was advanced in order to describe the low-energy properties 
of the 1-D HM [13]. In spite of the importance of these developments, a general description 
of spectral properties at finite frequencies is still lacking. 

We present in this Letter spectra for one-particle, spin- and charge-density excitations ob- 
tained with the maximum entropy method (MEM) [14] from quantum Monte Carlo (QMC) 
simulations of the 1-D HM. For the first time finite-frequency spectra are presented for sys- 
tem sizes up to 84 sites, allowing, irrespective of doping, for a clear identification of finite-size 
effects in different quantities . For half-filling, the one-particle spectra show besides the insu- 
lating gap, bands that closely follow a cosine-like dispersion. The spectral weight, however, 
is not evenly spread between the two bands, in contrast to mean-field calculations in the an- 
tiferromagnetic state. The magnetic structure form factor shows gapless excitations whereas 
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charge- density excitations show a gap. When the system is doped, a depression appears in 
the density of states at the Fermi-energy. A single band is observed, that follows very closely 
a cosine dispersion. In the doped case both spin- and charge-density excitations are gapless, 
with a "holon" -velocity larger than the "spinon" one, such that charge-spin separation is 
manifest. The spinon and holon velocities obtained lead to a value of the conformal charge 
of c = 0.98 ± 0.05 for N = 84 sites, in very good agreement with the exact value c = 1 
[12]. Finally, it is shown, that the spectrum for magnetic excitations fulfills excellently a 
frequency sum rule [15] for each &-point except for those close to 2kp, where departures of 
at most 10% are observed. 

We consider the 1-D HM described by the following Hamiltonian: 

H = -tJ2 (cj+i,*^,* + h.c?j + U J2 n i\ n il , (0-1) 

«,<r i 

where c^l are annihilation (creation) operators for an electron at site i with spin a, and 
m a = c\ a Ci )(J . Systems with a number of sites N ranging from N = 12 to N = 84 with 
periodic boundary conditions were simulated for inverse temperatures f3 = 1/ksT up to 
f3 = 20/ 1 and interaction strength U = At. The simulations were performed with the grand 
canonical algorithm [16], where the smallest values of At, i.e. the time-slice, was between 0.1 
and 0.125. The analytic continuation of the data to real frequencies was performed with the 
MEM [14], where the only "prior knowledge" used was the positivity of the spectral functions. 
We have chosen an uninformative default model ra 4 - = e, where e is a small quantitiy that 
merely suppresses noise in regions of insufficient information. The regularization parameter 
has been determined selfconsistently by classical MEM. For further technical details we refer 
to previous applications of the MEM to the single impurity Anderson model [17], the one- 
[18] and two- [19] dimensional Heisenberg antiferromagnet. 

We first compare the density of states D(u) of a ring with 12 sites between data from 
QMC simulations and exact diagonalization (ED), both for half-filling (Fig. la) and for the 
doped case (Fig. lb) with < n >= 0.833 ~ 5/6. Although the MEM is not able to resolve 
the rich structure obtained in ED, there is a good agreement in the shape and distribution of 
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weight in the half-filled case. In the doped case, the MEM can reproduce well the structures 
close to the Fermi-energy but has the general tendency of shifting weight to lower energies 
as one goes to higher energies due to the Laplace-transformation kernel. In particular, the 
pseudogap observed in ED that separates those states, whose weight was transferred upon 
doping from the upper Hubbard band (UHB) to the top of the lower Hubbard band (LHB) 
[20], from the remaining of the UHB, does not appear in the MEM data. 

In the following the data obtained for our largest system (N = 84), well beyond the 
capability of ED are shown [21]. We discuss first the spectral properties at half-filling. Figure 
2a shows D(u) and the spectral function A(k,u) is shown in Fig. 2b, where a dispersion 
can be clearly seen. Figure 2c shows the location of the maxima of A(k,u) with the errors 
assigned by the MEM. These errors give the uncertainty in the location of those maxima. The 
full curve is just a cosine function with a band-width that fits the QMC data. Remarkably, 
the dispersion obtained closely follows the two bands that would be obtained in a mean-field 
type of calculation, where antiferromagnetic order is assumed. However, the spectral weight 
is shared quite differently, since although most of the weight of the LHB appears for k<^, 
and the weight of the UHB is mostly concentrated in the region k <; -|, around k = -| weight 
is splitted between the two bands. 

Figure 3 shows the imaginary part of the spin-susceptibility xs(k } io) and the dispersion 
extracted from the maxima of Fig. 3a. From the slope of the dispersion around k = 0, 
we calculate a spinon- velocity v s /t = 1.23 ± 0.11. The exact value extracted from BA [7] 
vf A = 2tli(2ir / U) / Io(2ir / U) , where I and I-y are modified Bessel functions, is for U/t = 4, 
vf A /t = 1.2263 [22], that compares remarkably well with our value. The full line in Fig. 3b 
corresponds to the dispersion extracted from Bethe Ansatz [7] showing that the agreement 
with the exact results in the thermodynamic limit is not limited to the low-energy limit 
but extends to finite frequencies. For reasons of space, the equivalent quantities to those in 
Fig. 3 but now for charge-density excitations, are not displayed. We observe a rather broad 
structure with most of the weight around U, with a gap for A; — > that agrees within the 
errors with the one obtained in the one-particle spectrum. 
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Now we consider the doped case, where for definitness we have chosen a hole doping 8 ~ |. 



Figure 4 shows D(u), A(k } io) } and the dispersion of the maxima in A(k,u). A depression in 
D(u) is observed at the Fermi-energy ep, that we assign to the fact that a Luttinger liquid 
has a density of states D(u) ~| u |", with a ~ 0.038 for U/t = 4 [23]. The vanishing of 
D(u) at the Fermi energy cannot be seen in the simulation, since due to the small power a, 
an unrealistically high resolution would be necessary. Again, the "band-structure" obtained 
from A(k,u) can be very well described by a cosine band within the error bars. It should 
be stressed, that A(k,u) shows a peak that sharpens as one approaches ep for all the sizes 
simulated (N = 12,22,26,36,50,60,70 and 84). This seems to contradict the common lore 
of a Luttinger liquid having vanishing spectral weight at ep. However, Sorella and Parola 
have shown [11] in the limit U — > oo and for low density, that the quasiparticle weight 
z(kp) vanishes as A^ -1 / 8 , where the power coincides with a, the exponent characterizing the 
momentum distribution function around kp and D(u) around ep [23]. For the one-hole case, 
on the other hand, they obtained z(kp) ~ N -1 ^ 2 . Therefore, we should expect for U = 4 an 
even lower exponent [24], and hence, an utopically large system is needed to decide whether 
the system is a Fermi-liquid or not, on the basis of A(k,u). At U=4 we still do not observe 
significant weight in A(k,u) due to the UHB. However as was discussed in connection with 
ED, this can be due to the shifting of weight to lower frequencies by the MEM. 

The spectra for spin-excitations are shown in Fig. 5. Again, a linear dispersion appears 
for low energies around k = (inset of Fig. 5). Unfortunately the same resolution is not 
obtained around 2kp = bir /6. However, these data can be combined with the corresponding 
ones for charge-density excitations (Fig. 6). With the values of spinon and holon velocities 
extracted from the dispersion curves, it is possible to obtain the conformal charge of the HM, 
that is known exactly to be c = 1. Such a relationship stems from the conformal invariance 
of the model and has the following form [12]: 



where E (N) is the ground-state energy of the system with A^ sites, e is the ground-state 
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energy per site in the thermodynamic limit, v s and v c are the spinon and holon velocities, 
respectively, and c is the conformal charge of the model. The value obtained for N = 84 is 
c = 0.98 ± 0.05, where E and e were determined from Bethe-Ansatz [4]. Such a stringent 
test shows that QMC together with the MEM are able to give reliable real frequency data for 
low-lying excitations of the HM. It should be stressed, that the insets in Figs. 5 and 6 make 
manifest charge-spin separation in the model, since v c > v s . This is clearly seen in all system 
sizes simulated. Finally, we consider an exact sum-rule for the first frequency-moment of 
Imxs(&,cu) [15]: 

— co lm X s(k,co) = -< [[H,S'(k)],S'(-k)] > 

-oo Z7T 

= -2t(l - cosk) < H kin > , (0.3) 

where S z = |(^| ~ n i) an d < Hkin > is the expectation value of the kinetic energy, that can 
be accurately calculated by QMC. We find [21] that the sum rule is fulfilled accurately over 
most of the Brillouin-zone (deviations of less than 1%), with the exception of k > 2kp, where 
deviations (~ 10%) are obtained. They are probably due to a broad continuum similar to 
the one present in the 1-D Heisenberg antiferromagnet around k = tt [13,26]. This result 
together with the conformal charge obtained, demonstrate the degree of reliability of the 
numerical data from low to intermediate frequencies and for very large systems that are well 
beyond the capability of other methods like ED. 

Summarizing, we have presented real-frequency spectra for one-particle, spin- and charge- 
density excitations in the 1-D Hubbard model both in the insulating and in the metallic 
phases. The one-particle excitations show a band-like dispersion that can be accurately 
fitted by a cosine function. At half-filling, spin-excitations are gapless, whereas charge 
excitations show in this case a gap. The dispersions obtained agree with results from BA 
within the error bars. In the doped case, charge-spin separation is obtained, where the spinon 
and holon velocities excellently agree with the exact value for the conformal charge of the 
Hubbard model. The quality of the spectral data at intermediate frequencies is checked by 
an exact sum-rule, showing that only around 2kp, a departure of around 10% is present. 
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FIGURES 

FIG. 1. A comparison between QMC/MEM (lines with dots) and ED (full line) of the density 
of states for N = 12 sites: (a) Half-filling; (b) < n >= 5/6 (fi = 16, U = 4). 

FIG. 2. One-particle excitation for 84 sites at half-filling: (a) Density of states; (b) A(k,uj); (c) 
maxima of A(k,u), the solid line is a cosine function adjusted to the bandwith. (fi = 20, U = 4) 

FIG. 3. Spin-excitations at half-filling: (a) Imxs(k,u); (b) maxima of Imxs(k, u), the dotted 
line gives the slope at k = (v s /t = 1.23 ± 0.11), and the full line stems from Bethe-Ansatz. 

FIG. 4. The same as Fig. 2, but now for < n >= 5/6. 

FIG. 5. Imxs(k,tjj) for < n >= 5/6. The inset shows the linear part of the dispersion around 
k=0. The slope gives v s /t = 1.33 ± 0.09. 



FIG. 6. The same as Fig. 5 but for the charge susceptibility (v c /t = 2.02 ± 0.14). 
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